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Abstract: Functional renormalization group equations are analytically continued from imaginary 
Matsubara frequencies to the real frequency axis. On the example of a scalar field with 0{N) 
symmetry we discuss the analytic structure of the flowing action and show how it is possible to derive 
and solve flow equations for real-time properties such as propagator residues and particle decay 
widths. The formalism conserves space-time symmetries such as Lorentz or Galilei invariance and 
allows for improved, self-consistent approximations in terms of derivative expansions in Minkowski 
space. 
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1 Introduction 

In recent years, the functional formulation of Wilsons renormalization group [1] has developed from 
exact but rather formal relations such as the flow equations of Wegner and Houghton [2], Polchinski 
[3] or Wetterich [4] into a powerful tool for practical applications. By allowing for unconventional 
approximations and expansion schemes this made many interesting and partly non-perturbative 
investigations of quantum and statistical field theories possible. The field of applications ranges 
from cold atomic gases via QCD to Quantum Gravity, for reviews see [5]. 

In the modern formulation due to Wetterich one follows how a variant of the quantum effective 
action, the generating functional of the one-particle irreducible Feynman diagrams, changes while 
an infrared regulator is continuously removed. This evolution is described by an exact low equation 
of one-loop form. For a suitable chosen cutoff function only fluctuations of modes at a particular 
momentum scale k contribute to the flow. This organization of the field theory with respect to 
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the momentum scale has many advantages both for the intuitive interpretation and for practical 
calculations. 

So far, the formalism has been mainly applied in Euclidean type space or space-time, i. e. either 
static, classical statistical field theories where the fields depend on spatial position only or quantum 
field theories in the Matsubara formalism where the time and frequencies are imaginary. In all these 
cases one has the important practical advantage that propagators and other correlation functions 
have usually at most a single isolated singularity at vanishing frequency and momentum. (The 
situation is more complicated for fermionic theories at non-zero density where singularities appear 
for example on the Fermi surface.) 

In principle one can easily extend the formalism to Minkowski type space-times. Most of 
the definitions and derivations get modified only in a trivial way. However, in praxis a number 
of problems appear. Given that the invariant combinations —p^ + in the relativistic case and 
-po+f/{2M) in the non-relativistic case are not positive definite, the question arises which modes 
actually correspond to the infrared and which to the ultraviolet. It is therefore not straight-forward 
to construct appropriate regulator functions, in particular if space-time symmetries such as Lorentz 
or Galilean boost invariance should be preserved. Also, while for Euclidean type spaces it is well 
understood that the flowing action approaches the microscopic or classical action for large cutoff 
scales, this is in general not necessarily the case in Minkowski type space-time. 

Despite these problems, renormalization group methods have been applied to dynamical prob- 
lems. One example are classical, stochastic, react ion- diffusion problems [6]. This is facilitated by 
the fact that the microscopic inverse propagator is of the diffusion type dt — DV^ or iu + Dp^ in 
momentum space. Thus the propagator has no pole except ioi u = = and one can work for 
example with a regulator that depends on spatial momentum, only. 

A functional renormalization group formalism on the Schwinger-Keldysh closed time contour 
has been developed in refs. [7] and [8]. The goal of [8] was the investigation of non-thermal fixed 
points. Many technical problems - for example specifying a proper regulator function in Minkowski 
space - could be avoided in [8] by integrating the fiow equation and using a expansion to 
determine scaling relations between different correlation functions at a non-equilibrium fixed point. 
A different approach has been followed in ref. [7] to access far-from-equilibrium quantum field dy- 
namics. The regulator function was constructed with respect to time such that the renormalization 
group evolution and the time evolution coincide. 

Recently, functional renormalization has been used to investigate quantum gravity with Lorentz- 
ian signature [9]. The infrared regulator was chosen to depend on spatial derivatives only while 
the fiuctuations in the time directions have been regulated by compactification. A large amount 
of literature exists also on applications of different renormalization group techniques to transport 
through a zero- dimensional quantum system such as a quantum dot which is coupled to appropriate 
reservoirs [10]. 

Finally, real-time properties have been calculated from functional renormalization by analytic 
continuation. In refs. [11] and [12], the truncated renormalization group equations for a two- 
dimensional gas of non-relativistic bosons have been solved in the Matsubara formalism with imag- 
inary frequencies. The result for the propagator at the macroscopic scale k = has then been 
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analytically continued to real frequencies using numerical Pade approximant techniques. This al- 
lowed to calculate for example spectral densities and quasi-particle decay widths. A very similar 
approach has been followed in ref. [13] to calculate spectral functions of an impurity in a Fermi sea 
undergoing a transition from a polaron to a molecule. 

Note that such a procedure needs generically a rather large numerical effort. The reason is that 
the Pade approximant technique works only if enough (numerical) information about the dependence 
of the Propagator on the imaginary Matsubara frequency is provided. In other words, one has 
to calculate by some means the full frequency dependence of the (Matsubara) propagator. The 
numerical effort is even larger in a self-consistent scheme where the complete frequency-dependent 
propagator is used on the right hand side of flow equations. Nevertheless, such numerical approaches 
have been successfully implemented in the past, see for example ref. [14] for a scheme that allows to 
resolve both the momentum and field dependence of the propagator. This scheme has been applied 
to relativistic scalar fields at non-zero temperature in ref. [15]. 

An alternative is a non-self-consistent scheme where a simplified form of the propagator - ap- 
proximated for example by a derivative expansion - is used on the right hand side of flow equations. 
In that case the interesting and non-trivial information about the spectral properties are only ex- 
tracted from the flow equation but do not help to improve the results for other observables, e. g. 
thermodynamic potentials. 

In the present paper we develop an approach to the calculation of dynamic properties from func- 
tional renormalization which deviates substantially from all others that have been followed before. 
We use a linear response framework where the analytic continuation from imaginary Matsubara 
frequencies to real frequencies is done on the level of the flow equations and not on the level of the 
final result at the macroscopic scale k = 0. This has a number of advantages: 

(i) Since the flow equations of objects such as the propagator are usually available in analytic 
form or as an integral expression one can do the analytic continuation by hand and does not have 
to use involved numerical techniques such as Pade approximants. 

(ii) The interesting and non-trivial information about the real-time properties such as quasi- 
particle decay widths can be used in a self-consistent way on the right hand side of flow equations. 
This allows to improve the performance of a truncation in general so that also properties not 
directly connected to the propagator - for example thermodynamic properties - can be calculated 
with higher accuracy. 

(iii) Our formalism conserves translational symmetries as well as space-time symmetries such 
as Lorentz or Galilei boost invariance. Due to a convenient choice of the infrared regulator function 
it is nevertheless possible to perform the Matsubara summations in loop expressions analytically. 
This leads to well behaved expressions on the right hand side of flow equations where at most an 
integral over spatial momenta remains to be done numerically. 

(iv) In comparison with approaches based on the Schwinger-Keldysh closed time contour, the 
approach taken here is less involved. Also, since it is based on the well understood formalism in 
Euclidean space, one can profit from available knowledge. For example, it is known how the flowing 
action approaches the microscopic action for large cutoff scales or how useful regulator functions 
can be constructed. 
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On the other side it is clear that a formahsm that is based on analytic continuation is restricted 
to close-to-equilibrium situations. More specific, one can only calculate what can be accessed from 
linear response theory. More complicated non-linear response properties as they dominate in far- 
from-equilibrium situations are beyond the scope of this setup. 

This paper is organized as follows. In section 2 we recall the definition and the analytic prop- 
erties of the quantum effective action, the generating functional of the one-particle irreducible 
Feynman diagrams. This is followed by a discussion of the main principles of functional renormal- 
ization in its modern formulation for Euclidean space in section 3. Our main results concerning 
the analytic continuation of fiow equations are presented in section 4. We first discuss the principle 
idea and propose a suitable class of infrared regulator functions. On the example of a scalar field 
with 0{N) symmetry we explain then the technical details and show that the formalism works well 
in praxis. We solve the fiow equation within a truncation of the space of possible action functionals 
putting particular emphasis on real-time properties such as propagator residues and decay widths. 
Finally, we draw some conclusions in section 5. Appendix A collects some definitions of integral 
functions that are used throughout section 4. 

2 The quantum effective action and its analytic structure 

In this brief introductory section we recall the definition of the quantum effective action or one- 
particle irreducible effective action and discuss its properties, in particular its analytic structure. 
This section does not contain new results but combines existing knowledge that is usually discussed 
in different contexts. We use it to introduce our notation and to provide the basement for the 
discussion in later sections. 

2.1 The Schwinger functional 

The formalism we develop in the following can be applied to all variants of quantum field theories. 
In general the spectrum might contain bosonic and fermionic fields of different spin and in various 
representations of internal symmetries. However, in the present paper we are not interested in a 
particular model but in more conceptual issues. We concentrate therefore on the technically simplest 
case of a single scalar field. The presented formula can be generalized in a straight forward way to 
more complicated situations. (For fermionic Grassmann fields one needs some care in tracing the 
additional minus signs.) 

We are interested in the calculation of expectation values and correlation functions of the type 



Here, the operator 0[ip] may be a quite general functional of the field ip, for example a product of 
fields at different space-time points. The action 5* is taken to be real (hermitean) as appropriate 
for a quantum field theory after analytic continuation from Minkowski to Euclidean space. Usually, 
the action is given as an integral over a Lagrange density which depends on the fields (p and its 
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derivatives at the point x, 

S= I C{^{x),dMx),...). (2.2) 



To calculate expectation values and correlation functions as in Eq. (2.1), one introduces the partition 
functional 

Z[J] = [ Dip e-^^^^^^^^^. (2.3) 



By taking functional derivatives, one obtains expectation values 



as well as correlation functions 



We note that these objects depend on the source field J. 

Connected correlation functions can be obtained more direct from the Schwinger functional 

W[J]=\nZ[Jl (2.6) 

for example 

{ifixMy)). = {^{x)^{y)) - mm = ^^^t-^]- (2.7) 

For free theories, where S[ip\ is quadratic in the fields ip, one can usually calculate the Schwinger 
functional analytically. However, for non-trivial theories with interactions, this is not the case. If the 
interaction parameter is small, one can use a perturbative expansion. Various other approximative 
methods such as expansion or e-expansion exist, as well. 

It is crucial that W[J] encodes the information of all correlation functions. In terms of pertur- 
bation theory it contains contributions from all orders in the coupling constant as well as from all 
loop orders. This indicates that W[J] is usually hard to determine in a closed form. On the other 
side, a theory is basically "solved" if this can be done and many observables follow from functional 
derivatives of iy[J]. 

For practical purposes it is often useful to work not directly with the Schwinger functional W[J] 
but with a close relative, the quantum effective action T[(j)\. Its definition and properties will be 
reviewed in the following subsection. 

2.2 Quantum Effective action 

The quantum effective action is defined as the Legendre transform of the Schwinger functional 



r[0] = / J<P-W[Ji (2.8) 

J X 

with (p the expectation value 

m = j^,w[J] = (^(^))- (2-9) 

dJ[x) 
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Directly from its definition one obtains the field equation 



S 



m 



J{x). 



(2.10) 



Taking another derivative, one finds 



m 



{W^'^[J])-\x,y). 



(2.11) 



54){x)S(l){y) 



S(l>{x) 



The quantum effective action is useful for several reasons. Since it is the Legendre transform of 
the Schwinger functional, it contains basically the same information. This implies that correlation 
functions of all orders in perturbation theory can be obtained from functional derivatives of r[0]. 

In perturbation theory r[0] can also be defined as the generating functional of the one-particle 
irreducible Feynman diagrams. While this seems to be a quite formal statement, it has a direct 
practical consequence: The quantum effective action is exact at tree level. By this we mean that 
if the quantum effective action r[0] is known, everything that remains to be done to calculate a 
physical observable such as a scattering cross section is to evaluate tree-level Feynman diagrams with 
the vertices and propagators taken from the quantum effective action r[0]. This implies in particular 
that the quantum effective action contains directly the physical parameters such as masses, charges, 
magnetic moments and so on as they arc measured in experiments. No loop diagrams need to be 
calculated any more. The parameters appearing in r[0] are already the renormalized ones. It is 
important to emphasize at this point that the quantum effective action contains not only vertices 
that are "renormalizable" in the sense of perturbative field theoretic renormalization theory. In 
contrast, it contains in general all terms that are allowed by symmetries. For example, the quantum 
effective action of QED will contain a term of the form 



with a coefficient k that is related to the anomalous magnetic moment g — 2 hj k = ^{g — 2). 

To see that F [cp] is "exact at tree level" we note that the Schwinger functional can be obtained 
from the quantum effective action by a Legendre transform 



see in different ways that the exact connected n-point functions generated by W[J] can be written 
as tree diagrams with propagators and vertices obtained from F[0]. 

One possibility is a "proof by example". Taking functional derivatives of W[J] one finds (we 




(2.12) 




(2.13) 



where (J) on the right hand side is determined from the implicit equation J 



AF[0]. One can now 
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use a symbolic notation for simplicity) 

^(1) = —W[J] = K 

(2) _ J_ J_^r .p(2)wl 

" <5J« 5J^ ^ ~ 6J^ ~ (2.14) 

w(3) _ _/p(2)\-l fr{2)N-l ^p(2)wlp{3) 

w'Xs = ■■■ 

By continuing this one finds indeed that W^"^ can be written as a tree level expression involving 
the propagator (r*^^))"^ and the vertices T^^\ 

Another, more elegant way to see this [16, 17] is to consider for a moment a functional integral 
where r[(p] replaces the microscopic action S[ip] 

^Wr[J]^ j D^e--^^^^^+-^^'^. (2.15) 
A loop expansion for VFr[</] is an expansion in the parameter g with 

(W^r[J])^,oops~/~'- (2.16) 

In particular, the term of order g^^ contains only tree diagrams. One evaluates now the functional 
integral in Eq. (2.15) in the limit (? — ?■ 0. It is then dominated by the first term in the saddle point 
approximation, the stationary phase term 

\im{gWM)= I Jv-m. (2.17) 

with ip determined by J = 5V[(p\/5ip. However, this is just the defintion of W[J] in Eq. (2.13). 
In other words, we found that the Schwinger functional W[J] can be obtained from a "functional 
integral" as in (2.15) where the microscopic action S is replaced by F and only tree diagrams are 
taken into account. 

2.3 Perturbative loop expansion 

It is clear from the considerations in the previous section that the quantum effective action is a 
useful object. On the other side it is also clear that it is general not easy to calculate. To investigate 
the connection between r[0] and the microscopic action S[(p\ more closely, it is useful to consider 
the following implicit functional integral representantion which can be derived directly from the 
definitions after a shift in the integration measure 



e 

We have restores here the units of h for a moment. One can now understand the perturbative loop 
expansion as an expansion of the right hand side in Eq. (2.18) in h. One finds 

T[(P\ = const + S[^] + n^STrln5(2)[0] + . . . (2.19) 
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The ellipses stand for two loop and higher order loop expressions. Assuming that the terms mul- 
tiplying powers of h on the right hand side are finite, one finds T[(f)] = S[(l)] (up to an irrelevant 
constant) in the classical limit 0. This is the reason why S[(j)\ is sometimes called "the classical 
action" . 

2.4 Analytic structure 

So far we have restricted the discussion to Euclidean quantum field theory ("imaginary time"). 
Concepts like the Schwinger functional or the quantum effective action are of course also useful for 
accessing real time properties. We note at this point that the definition of the effective action T[(j)] 
can in principle be easily extended from Euclidean to Minkowski space, at least in the vacuum for 
vanishing temperature. The fields (p of the functional integral depend then on the real time variable 
t in addition to the space position x. The correlation functions one calculates as expectation values 
of field monomials correspond to time-ordered correlation functions in the operator picture. 

However, we take here another point of view where the configuration space for the functional 
integral fields (f and their expectation values is always Euclidean space. (At non-zero temperature 
it is a generalized torus with circumference 1/T in the imaginary time direction.) ffowever, the 
n-point functions which can be obtained from T[(f)] as functional derivatives, can be analytically 
continued towards the real frequency axis. For example, from the second functional derivative of 
r[0] with respect to the fields one can obtain the propagator G{po,p) at the imaginary Matsubara 
frequencies po — ioUn- As will be discussed below, it follows from general principles that this function 
can be analytically continued to the complete plane of complex frequencies except for possible poles 
and branch cut singularities on the real axis. 

A similar analytic continuation is possible for the higher order correlation functions and in 
this sense the effective action in Euclidean or Matsubara space contains already all the relevant 
information for the dynamics in Minkowski space. It is clear, however, that within the approach 
propagated here one can only access closc-to-cquilibrium properties in the spirit of linear response 
theory. In contrast, a more elaborate setup based for example on the Schwinger-Keldysh formalism 
would be needed in a far-from-equilibrium situation. 

Let us now discuss constraints on the analytic structure of n-point functions as they arise from 
basic principles of quantum field theory in more detail. For simplicity we restrict the considerations 
here to the two-point function F*^^) but we stress that at no point we assume that the fields wc 
consider are fundamental. Therefore, by applying similar arguments to two-point correlations of 
composite fields one can in principle cover much of the analytic structure of a general n-point 
function F^"). In a homogenous situation, the two-point function can be written as 

.f') = ^ = (2.)'f "(P - P') W (2.20) 

with G{p) being the momentum representation of the Euclidean propagator. 
For a microscopic action S[if] of the standard form 

S[ip] = J d^x C{ip,d^ip), (2.21) 
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one can go from the functional integral representation to an operator picture. Translations with 
respect to (real) time are generated by a Hamiltonian with real and positive eigenvalues. Similarly, 
momentum operators with real eigenvalues generate spatial translations. The Hilbert space gets 
spanned by a complete set of states which can be taken to be eigenstates of the four-momentum 
operator P^. Using this formalism as well as requirements from space-time symmetry, causality and 
unitarity, one can derive the following Kallen-Lehmann spectral representation of the propagator 
(for a detailed derivation see e.g. [17]) 

/•OO 1 

G{p)= / rf^^p(/i^)^— ^ (2.22) 

with real and non-negative spectral weight p(/i^) > 0. 

The representation (2.22) is valid for an Euclidean quantum field theory where = jp' > 
but can also be used for the analytic continuation to Minkowski signature. In Minkowski space 
one has to pay attention to the singularities at Pq = + /i^ and specify an integration contour in 
the complex frequency plane. Often this is done by introducing additional infinitesimal terms ie to 
move the singularities slightly away from the real po-axis along which one integrates by convention. 
In the representation 

dfi^ p(/i2) ^ I ^ , ^ ) (2.23) 

+ ii'^ y-po + ^/jP + fi'^ ±ie -po - ^/ fp + fi"^ ± ie J 

the combination —ie, +ie in the first, respectively second term corresponds to the time-ordered 
propagator while the combinations —ie, —ie represents the retarded, +ie, +ie the advanced and 
+ie, —ie the anti-time-ordered propagator. However, this ie-prescription is just one way of doing 
the book keeping and deforming slightly the contour itself is just as good. 

We emphasize again that the derivation of the spectral representation (2.22) works only for 
microscopic actions of the standard form (2.21). It is important that it is local and that the 
Lagrangian depends only on the fields as well as its first derivatives d^ip. In particular, if one 
adds an infrared regulator term of the form 

A^.M = J ^,lv{-p)Rk{pMp) (2.24) 

to the microscopic action, this corresponds in general to a non-local action, or a situation where 
arbitrary high derivative orders appear in the Lagrangian. It is not clear how an operator formalism 
can be developed in this case. One can choose the regulator function such that it respects 
space-time symmetry and causality requirements but, nevertheless, one can in general not expect 
that the propagator corresponding to a theory regularized in this way admits a standard spectral 
representation. 

A few important properties can be read off from Eq. (2.22). First, for Euclidean signature or 
p^ > the propagator is real and positive. This reflects a general property of Euclidean quantum 
effective actions closely related to Osterwalder-Schrader reflection positivity. In contrast, close to 
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the real frequency axis where Im po = one can write {V stands for principal value) 



G{p) = / dfj,^ p{fJ'^)'P — 2,^1 — 2 + ^ ^ sign(Re po) sign(Im po) p{pl - f). (2.25) 

This expression makes explicit that there is a branch cut on the real p^ axis where the imaginary part 
of G{p) changes its sign for those combinations of frequency and momentum where p{p'^ — fp') > 0. 
Away from the real frequency axis the propagator G{p) is analytic and, since p(/U^) cannot vanish 
for all /x^, also non-zero. 

In summary, the propagator G{p) as a function of the complex frequency po has poles, zero- 
crossings and branch-cut singularities only on the real axis. The inverse propagator G{p)^^ has 
according to the above discussion all its zero-crossings, poles and branch-cut singularities on the 
real frequency well. 

Moreover, from the representation in Eq. (2.25) one can see that close to the real frequency 
axis one has 

P{p) = G{pr' = P,{pI - f) - ts{po)P2{pl - f) (2.26) 
where we have introduced the abbreviations 

s{pq) = sign(Re po) sign(Im po), (2.27) 

P^iPl-p") = rf 2 ^,12 (2-28) 

[Re G(p)]2 + [vr p{pt, - 

and P2{Pq — P^) > is related to p{j)Q — fP') by 



^^(^^0 - f) = rp r( ill f 2 ^^12 • (2-29) 
\ReG{p)Y + ['K p{pi-p^)Y 



Close to a point Pq — p^ = fnP where Pi vanishes, one can use an expansion of the form 

Pi = Z(-p2+pe + m2) + ... 
P2 = + . . . 

such that for 7^ = mV and T <^m the propagator is of the Breit-Wigner form 



(2.30) 



Z {-pI + p' + m^y + m^T^ " ^ ^ 

The representation (2.26) and expansions of the form (2.30) will be widely used in this work; 
not only for the full inverse propagator P{p) but also for a scale-dependent generalization thereof. 
Before discussing that in more detail in section 4, we recall some important concepts and principles 
of the functional formulation of the renormalization group for Euclidean space in the next section. 
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3 Principles of functional RG in Euclidean space 



In this section we recall the main ideas and principles of functional renormalization for Euclidean 
quantum field theories. The flowing action or effective average action is introduced as a general- 
ization of the quantum effective action. We discuss how it interpolates continuously between the 
microscopic or classical action and the quantum effective action while an infrared cutoff is removed. 
This evolution is described by an exact functional differential equation. Equally important for prac- 
tical applications are approximative methods to solve it and we discuss the principle ideas behind 
these approximations towards the end of the section. 

The present section is kept rather brief and introduces only the prerequisites for a later dis- 
cussion of analytically continued flow equations in section 4. For a more detailed discussion of 
functional renormalization in Euclidean space and applications to many concrete problems we refer 
to many good review articles in the literature [5]. 

3.1 Exact flow equation 

We start by modifying the partition functional and Schwinger functional introduced in (2.3) and 
(2.6), respectively, according to 

where we have added an infrared cutoff term 

ASfeM = \jv{-p)Rk{p)v{p) (3.2) 

to the microscopic action S[ip\. The function Rk{p) is real and positive for Euclidean momentum 
argument p. It can be of different form in different situations and we discuss possible choices below. 
As a general feature, one requires that Rkip) serves as an infrared cutoff as the scale k which is 
usually implemented by the property 

lira Rk(p) ^ k"^. (3.3) 

p-s>0 

In addition, Rk{p) should suppress fluctuations on all momentum scales for very large values of the 
cutoff scale k. This requires Rk{p) to diverge for — oo. Usually this divergence is taken is taken 
as strong as k"^. Finally, the cutoff function should vanish for /c — > for all p, 

limi?fc(p) = 0. (3.4) 

It is now convenient to define a modified form of the effective action. The flowing action or 
effective average action F^, [0] is defined by subtracting from the Legendre transform 



ffc[0] = / J<P-Wk[J] 



(3.5) 



with = j^l4^fc[J] the cutoff term 

r,[0] = ffc[0]-A5feM. (3.6) 

Using the above definitions, one can obtain an exact functional renormahzation group equation 
for the flowing action rfc[0], the Wetterich equation [4] 

duTk[<p\ = ^Tr(rf [0] + RuY^duRk- (3.7) 

We note that the right hand side has the form of a one-loop expression. Nevertheless, it is an exact 
equation taking all orders of perturbation theory as well as non-perturbative effects into account. 
Other properties of [</>] follow from an implicit functional integral representation 

g-r^M = j Dip e-'5'[<^+'^]-^-^feM+/(s|;rfc[</'])v'_ ^3_8) 

Similar to (2.18) this equation follows directly from the definitions together with a shift in the 
integration measure. 

For k oo and therefore i?^ — )■ cxo one finds that the cutoff term ASk dominates the functional 
integral over the fluctuating field ip in (3.8). It leads to a strong suppression for large values of 
y^{—p)^{p) = f*{p)y^ip) and a saddle point expansion becomes a good approximation. One has 
therefore for very large k = A 

Ta[^] = const + S[^] + ^Trln(5(2)[0] + Ra). (3.9) 

For many applications the constant ((/)-independent) part is not important. 

In addition, the contribution from the one-loop term (the second term on the right hand side 
of (3.9)) is suppressed by negative powers of -Ra for derivatives of (3.9) with respect to the field 0. 
This implies that the flowing action effectively approaches the microscopic action for large values 
of the cutoff parameter k, 

lim r40] = 5[0]. (3.10) 

fc— >oo 

To summarize, the flowing action Tk[4>] has been defined as a modified Legendre transform of the 
Schwinger functional ^^^[J] in presence of the infrared cutoff function Rk{p)- For very large values 
of k it approaches the microscopic action S[(j)] (up to an irrelevant constant as well as possibly 
one-loop terms) and for = it equals the quantum effective action r[0]. In other words, the 
flowing action flows from the microscopic action S[(j)] to the quantum effective action r[0]. 

The exact functional differential equation for the flowing action (3.7) together with equation 
(3.9) can be seen as a differential formulation of the functional integral formalism. As an alterna- 
tive formulation of quantum field theory this has an advantage for some formal questions but the 
most important application is for concrete practical and often approximate calculations. We will 
briefly discuss such an application for Euclidean field theory below. The emphasis is on the main 
principles of the approximation scheme rather than on the numerical results. It is useful to have 
these principles for Euclidean space in mind when the analytic continuation to Minkowski space is 
discussed thereafter. 
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3.2 Derivative expansion 

The example we discuss here is a classical scalar field theory with 0{N) symmetry in d spatial 
dimensions (usually one has d = 3 but also d = 2 or d = 1 and even a single point with d = are 
of experimental relevance). The microscopic action S is of the form 

S = j d^X |^V0„V0n + ^mVn0n + (0n0n)'| ■ (3.11) 

This model can be used for example to describe the critical properties of many second order phase 
transitions. 

The cutoff term for the investigation of the model in Eq. (3.11) has the form (with 0*(p) 

(2-kY 



:-p1 and I = / ^; 



The function Rk{p) depends only on the magnitude p = so that the cutoff term is invariant 
under rotational symmetry. Note that it respects also other symmetries as the internal 0{N) 
symmetry and translation. Useful cutoff functions are for example the exponential cutoff 

Mp) = JlS (3-13) 

or the one proposed by Litim [18] 

Rk[p) = Zk(e-p')Q{e-v'). (3.14) 

Here is the wavefunction renormalization constant and one has Za = 1. As a feature of both 
functions in Eqs. (3.13) and (3.14) we note R^ij) = 0) = Zl^ and Rk{p) — ^ for p^/A;^ — )■ oo. One 
can see these functions as a p-dependent mass term for the modes with p^ ^k'^. 

That Rk{p) goes to zero for large values of has a number of practical advantages. For a 
fixed value of k, the fact that Rk{p) ~ for ^ fc^ means that fiuctuations of modes with large 
momentum are integrated in the flowing action Tk, already. If k is lowered further, there will be no 
change of the correlation functions due to fluctuations of these modes. In praxis, this means that 
the momentum integrals that appear on the right hand side of the flow equation of an object such 
as the effective potential have no (or only small) contribution from momenta with p^ ^ k'^. There 
is an effective UV cutoff due to the term dtRkip) on the right hand side of the flow equation. 

This picture can also be conflrmed from the functional integral expression for the flowing action 

e-rfeM = j DyD g-5[0+<^]-/^|<{p)/?,(p)^„(p)+/pl^^"(p). (3.15) 

For a given value of k, the modes of the fluctuation fleld ^Pn{p) with p^ ^ k"^ are not affected by 
the cutoff term. In contrast, for the modes with p^ ^ k"^ the cutoff term leads to an additional 
Gaussian suppression similar to a mass gap of size k"^. 
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With a cutoff function as in Eq. (3.13) or (3.14), the contribution of long-ranging modes with 
small values of are included in a steady way in the functional integral in Eq. (3.15). These are 
precisely the modes that become important close to second order phase transitions and dominate 
the universal critical dynamics. A cutoff function that decays fast enough for large values of /k'^ 
allows to separate the treatment of these modes from the higher momentum modes which are not 
relevant for the critical phenomena. 

The above described separation of cutoff scales reahzed by a cutoff function Rk{p) that decays 
sufficiently fast for large p makes an expansion in terms of spatial momenta a useful approximation 
scheme for many purposes. In position space this corresponds to a derivative expansion. One writes 

r,[0] = j d'^x |f/fc(p) + ^Zfc(p)V0aV0a + ^n(p)VpVp + 0(9^)| , (3.16) 

with p = ^(j)a4>a- The lowest level only includes the scalar potential and a standard kinetic term. 
The first correction includes the p-dependent wave function renormalizations Zk{p) and Yk{p). The 
next level involves then invariants with four derivatives and so on. 

On first sight there is no reason why a derivative expansion as in Eq. (3.16) should be a good 
approximation. In principle both small and large momenta contribute on the right hand side of 
the flow equation and even observables at small momenta could be affected by the large momenta 
regime where (3.16) is expected to converge badly. To understand why it nevertheless works rather 
well in praxis consider the following generalization of the scheme 

rfc[0]= I t/^x|f/fc(p) + ^V0a^fc(p,-V')V0a 

. (3.17) 

+ -vpn(p,-v2)vp + ...|. 

The functions and depend now also on the momentum p = \f^. The scheme in (3.17) has 
the advantage that it allows to resolve the full momentum dependence of the propagator. However, 
in praxis one usually neglects terms higher than quadratic in the momenta. Nevertheless, derivative 
expansion often leads to quite good results. The reason is the following. On the right hand side 
of the flow equation the cutoff insertion i?fc(p) in the propagator (r*^^)(p) + Rk{p))~^ suppresses the 
contribution of the modes with small momenta. On the other side, the cutoff derivative dkRk{p) 
suppresses the contribution of very large momenta provided that Rk{p) falls of sufficiently fast for 
large p. Effectively mainly modes with momenta of the order k"^ contribute. 

On the other side, if one evaluates flow equations for objects as Zk{p,f^), Yk{p,p^) etc., one 
main effect of the external momentum is to provide an infrared cutoff scale of order p"^. Such an 
infrared cutoff scale is already provided by Rk itself and one might therefore also work with the 
fc-dependent but p^-independent couplings 

Zk{p,p' = 0), Yk{p,p' = 0). (3.18) 

This brings us back to (3.16). In other words, a properly choosen cutoff function improves a 
truncation with ^-independent coefficients Zk and Yk such that some of the momentum dependence 



- 14 - 



is effectively taken into account. We emphasize that it is important that the cutoff Rk{p) falls off 
sufficiently fast for large p. If this is not the case, the derivative expansion might lead to erroneous 
results since the kinetic coefficients as appropriate for small momenta and frequencies are then also 
used for large momenta and frequencies. Only when the scale derivative dkRk{p) provides for a 
sufficient ultraviolet cutoff does the derivative expansion work properly. 



4 Analytic continuation of flow equations 

In this section we discuss flow equations in Minkowski space. The original formulation of the 
formalism as outlined in section 3 has been in Euclidean space or for imaginary times. There are 
now different conceivable strategies to obtain results about real-time propagators which will be 
discussed below. 

Our general point of view is that the flowing action Tk[(j)\ is a priori defined in the Matsubara 
formalism, i. e. within a functional integral setup where the configuration space for the fields (p 
and their expectation value has a cyclic imaginary time direction with circumference 1/T. The 
correlation functions that follow from Tk[4>] as functional derivatives can then be analytically con- 
tinued from the discrete, imaginary Matsubara frequencies to the complex frequency plane and in 
particular towards the real frequency axis. 

For some properties of the flowing action the analytic continuation has no effect. This is the 
case in particular for the effective potential which is obtained from the flowing action by evaluating 
it for fields that are constant in space and time so that analytic continuation is trivial. On the other 
side, for frequency- and momentum-dependent objects such as the propagator or other correlation 
functions, it does make a difference. They can either be evaluated in Euclidean space with imaginary 
frequency argument or - after analytic continuation to real frequencies - in Minkowski space. 

In most situations one is finally interested in macroscopic quantities for k = 0. One way to 
get for example the macroscopic propagator in Minkowski space would be to solve its flow equation 
with an arbitrary imaginary value go = "^^^ of the frequency argument (at finite temperature the 
frequency is restricted to the discrete values go = i'27cTn in the Matsubara formalism) and to do 
the analytic continuation only at the end, i.e. for k = 0. The advantage of this procedure is that 
the flow equation formalism is only used in Euclidean space where it is rather transparent and well 
understood. On the other side, the analytic continuation can be very difficult in praxis since the 
propagator in Euclidean space is usually available at best numerically. Methods based for example 
on Fade approximants need information from many Matsubara frequencies and the numerical effort 
gets rather large. Nevertheless, this approach has been successfully followed in the past and allowed 
to calculate for example the spectral function of non-relativistic bosons in two spatial dimensions 
[11, 12]. 

In this paper we propose another way to solve the problem, however. Instead of doing the 
analytic continuation only after solving the flow equation or at /c = 0, we will analytically continue 
the flow equations themselves. Since in contrast to their solution the flow equations are usually 
available in analytic form, it is not necessary to use involved numerical techniques. As an example 
we investigate in particular the the flow equation for the propagator. For particular choices of the 
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infrared cutoff function Rk we show how it is possible to evaluate this object for real frequency 
argument and to solve the flow equations directly in Minkowski space. No analytic continuation 
of the final result is then needed any more to access the real time properties. Moreover, it will 
become apparent that such a procedure also has some advantages in designing approximations. A 
derivative expansion in Minkowski space can be done as a Taylor expansion around the frequency 
corresponding to a pole or branch-cut singularity of the propagator. Loop expressions that appear 
on the right hand side of flow equations but also on-shell properties of the effective action itself are 
strongly dominated by these singular structures. One can therefore expect that the convergence 
properties of such an expansion in Minkowski space are better than those of a derivative expansion 
around vanishing values of the frequency as it can be done in Euclidean space. 

This section is organized as follows. In subsection 4.1 we discuss a class of regulator functions 
for which the analytic continuation of the flow equation to real frequencies can be done, at 
least for truncations based on derivative expansion. We apply this method subsequently to a model 
for scalar fields with 0{N) symmetry. The truncation for this particular system is introduced in 
subsection 4.2 while subsections 4.3 and 4.4 discuss the flow equations for the effective potential 
and propagator, respectively. Subsection 4.5 presents some numerical results for the solution of the 
flow equations within the truncation. 

4.1 Regulator 

In this subsection we discuss the particular choice of the infrared regulator function Rk that allows 
for analytic continuation of the regularized propagator (P^ + Rk)~^ from the imaginary frequency 
axis or Euclidean space to real frequencies. The discussion will be based on truncations where the 
inverse propagator is close to the real frequency axis of the form 



where Z, z, m? and 7^ are fc-dependent, real and positive quantities. 

Let us now choose a regulator function Rk{po,p)- It is a priori not clear what requirements must 
be fulfilled when choosing a cutoff function in Minkowski space. For example, it is not clear what 
the conditions for Rk{po,p) are to be a good infrared regulator nor is it obvious what is required 
for it to have good regulating properties in the ultraviolet. Actually, it is not even obvious which 
modes correspond to the IR or UV, respectively, given that the combination —pl+j? is not positive 
definite. 

On the other side, all these issues are well understood for Euclidean signature, see the discussion 
in section 3. The approach we follow here is therefore to choose a regulator function with all the 
desired properties in Euclidean space or for positive argument — + p^ > and to use analytic 
continuation to extend the calculations to Minkowski space. 

For many choices of Rk such a program would be rather difficult to perform since a function that 
is smooth and regular on the imaginary frequency axis may nevertheless have poles and discontinu- 
ities in other regions of the complex frequency plane. One can even argue that this is unavoidable 
if one wants Rk{po,p) to be a function that decays with large imaginary values of Pq. This in turn 




(4.1) 
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is of course necessary for Rk to serve as an effective UV regulator on tlie riglit hand side of the flow 
equation. 

To make the analytical continuation possible in praxis we choose a class of rather simple IR 
regulators with algebraic decay in the UV: 

Mpo,P) = Ze , , ^ , , • (4.2) 



fc2 j \ y 1,2 ^ 

The coefficient Zk can be chosen conveniently. Usually it is taken to agree with the wave function 
renormalization constant but more generally it can be some real and positive function of the scale 
parameter k with not too strong dependence on k. 

The function (4.2) has all the desired properties for Euclidean arguments —pi + > if the 
coefficients Cj are real and positive. It is clear that the UV regulating properties become better if 
some Cj with large j are non-zero. On the other side, practical calculations are simple if only a few Cj 
with small j are non-zero. Arguably the simplest, non-trivial choice is Ci = c > 0, C2 = C3 = • • • = 
and we will discuss it in more detail below ^. One should keep in mind, however, that the decay of 
Rk for large momenta and Matsubara frequencies is rather mild for this choice. For the convergence 
properties of the widely used derivative expansion it may be desirable to have a quicker decay, either 
with a higher power or exponential in the momentum argument. 

In the following we discuss the analytic structure of the regularized propagator [Pk + Rk)~^ as a 
function of a complex frequency argument Pq. In principle, this depends on the complete functional 
form of Pk and it is not sufficient to know the form of Pk for example close to the real frequency 
axis. On the other side, at least when k is sufficiently small one expects the relevant properties of 
the propagator to be dominated by the on-shell singularities, i. e. the poles and branch cuts on the 
real frequency axis. Following this rationale, we use the expression (4.1) not only close to the real 
frequency axis but everywhere in the complex p^ plane, keeping in mind that the coefficients Z, iv? 
and 7^ should be determined at the singularity, i. e. close to the point on the real frequency axis 
where the real part of Pk vanishes. 

With the above choice of the regulator function it is straightforward to show that the regularized 
propagator can be written as (we use fh? = m? /k'^ and 7^ = 7^/A;^) 

1 / /3i , P2 



with 



1 fl fn^ 72\ 
"1/2 = 9 ( - ^ -5(^0) — ±[A + i s{po) B] , 



Pi/2 = Y^±[C + is{po)D] 



(4.4) 



^ A suitable cutoff function for non-relativistic field theories is obtained by replacing — Pq with —pQ +p^/ (2M) 
in equation (4.2). 
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We choose the branch cut of the square root function to be along the negative real axis. The objects 
A, B, C and D are then always real. Moreover, using the above definitions it is straightforward to 
establish that for all values of c > 0, fh? > and 7^ > one has 



1 

2 \c 
B>0 



+ 



~ 1 
m 



± A > 0, 



(4.6) 



< 0. 



The decomposition in (4.3) is very useful for practical calculations since it resembles closely 
the form of a free propagator. A similar representation is actually possible for the whole class of 
regulators in (4.2). 

We note at this point that in contrast to the propagator P^^ in Eq. (4.1), the regularized 
propagator (P^ + Rk)~^ in (4.3) has singularities away from the real frequency axis. Indeed, for 
^ — B < and s{po) = 1 there are poles at po = ±^yfP + aikP'. A spectral representation as in Eq. 
(2.22) is therefore not possible for (P^ + Rk)~^- Although we do not give a proof here we believe 
that this is a generic feature for cutoff functions that serve as an effective UV regulator. 

One of the most important features of the regulator function proposed here is that one can 
calculate in praxis with the decomposition (4.3). Besides the pole singularities, (P^ + Rk)~^ has 
also a branch cut originating from the second term in (4.3). However, in an approximation where one 
assumes that all integrals along this branch cut are dominated by the nearby poles on the different 
Riemann sheets one can always perform the frequency integrals (or Matsubara summations for 
T > 0) analytically which is an enormous advantage for practical calculations. 

4.2 Truncation for 0{N) model of scalar fields 

In the remainder of this section we perform the analytical continuation for the flow equations of 
a concrete model, the 0{N) model for a scalar field in the phase with spontaneous symmetry 
breaking. The excitation spectrum consists of a massive radial mode as well as — 1 massless 
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Goldstone modes. Due to the non-linear coupling, the radial mode can actually decay into two 
Goldstone excitations which gives rise to a non-vanishing decay width. 

Since the breaking of the 0{N) symmetry is spontaneous, the effective action and its k- 
dependent generalization, the flowing action rfc[0] are invariant with respect to global 0{N) sym- 
metry transformations. It must also be invariant under translational and rotational symmetries as 
well as under Lorentz transformations. Apart from these constraints it may be a rather general 
functional of the field 0. However, to make a (numerical) solution of the flow equations possible it 
is in any case necessary to truncate the space of possible functionals rfc[0]. 

In this paper we use a truncation based on derivative expansion in Minkowski space. More 
speciflc, we take to be of the form 

t M=l 

with p = \ ^M=i (fu- some arbitrariness in this decomposition by demanding the momentum 

dependent parts to vanish for po = P = 0, i.e. 

P^(0,0) = P,(0,0) = 0. (4.8) 

We now make the crucial assumption that the momentum dependent propagator parts P(j,{po,p) 
and Pp{po,p), considered as a function of a complex frequency po, are of the same analytic structure 
as the full inverse propagator in (2.26). This is certainly the case at the macroscopic scale k = 
since rfc=o[0] = T[4>]- For non-zero k one expects in principle deviations due to the frequency 
dependence of the cutoff function Rk{po,p). One can see it as an element of our truncation that 
these deviations are neglected, however. 

In this spirit and taking constraints from Lorentz invariance into account we write the functions 
Pfj, and Pp as 

P^/pipo,p} = Z^/pi-pl+f) [-pl+f] - tsipo)f^/pi-pl+f). (4.9) 

We use again the function s{po) deflned in Eq. (2.27). The functions 7|/p(— Po + P^) nonzero 
only for negative argument. 

To derive the flow equation for tJk{p), we expand (4.7) around a background conflguration 
according to 

0i(t,x) =0 + 501 (t,x), ^Mit,x)=6^Mit,x) (M = 2,...,Ar) (4.10) 
and keep only terms that are quadratic in the "fluctuating parts" 50m- This yields (p = ^0^) 

Tk,2 = I ^50i(-Po, -P) [Pi>ipo,P) + pPp{po,P) + U'liP) + 2pU'^{p)\5(pi{pQ,p) 

N ^ _ _ _ _ (4.11) 

+ X] i^Hm{-Pq, -P) [P^{po,p) + Uk{p)]6(j)M{po,P) 

M=2 
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In principle one could now derive a closed set of flow equations for the set of functions Z^{—pQ + 0^), 
^pi~Po 7|(~Po 7p(~Po f/fc(p) by using appropriate projections of the flow 

equation for F^. At the macroscopic scale k = 0, the functions Z^{—pl + p^), Zp{—pl + p^), 
■y'^{—pl + f^) and 7p(— + V^) would contain all the interesting information of the propagator, 
spectral density etc. However, this procedure would result in a set of coupled, partial integro- 
differential equations which is rather hard to solve numerically. 

Instead we therefore devise here a simpler approximation based on the observation that the 
physical properties of a propagator are to a very large extend dominated by singularities, in par- 
ticular poles or branch cuts corresponding to zero-crossings of the inverse propagator or at least its 
real part. Around these singularities it is sensible to use a derivative expansion or Taylor expansion 
of the inverse propagator with respect to its frequency and momentum arguments. 

We observe at this point that the frequency po where the propagators become singular depend 
on the value of the background field p. The most important properties will be dominated by the 
minimum po of the effective potential where U'{po) = 0, however. We therefore choose for the 
inverse propagator of the Goldstone modes the expansion point po = p = 0. The imaginary part 7^ 
vanishes at this point and we write 

P^{Po,p)^Z^{O){-pI+0'), (4.12) 

the higher orders being neglected. 

The situation is different for the radial field. Since it is massive, it is sensible to choose the 
expansion point as po = mi and write there 

P<t>{Po, P) + PoPpiPo, V) ,^ 
^ [Z^{-ml) + -p,Zp{-m\)\ {-pI zs{po) [f^{-ml) + Pofpi-mD] , 

with higher order terms being neglected. 

We define now the wave function renormalization as Z^ = Z^{0), the renormalized fields as 
= \fZ~4>^M, P = Z^p, the effective potential Uk{p) = Uk{p) and the anomalous dimension as 

?70 = —-^dtZ^. In addition we introduce the abbreviations Zi and 7^ by 



(4.14) 



An expression analogous to (4.11) but in terms of renormalized quantities can now be written 
Tk,2 = I ^50i(-Po, -P) [Zi [-vl+f] - ts{po)ll + U'M + 2pUM] SMpo,P} 

•J P,PO 

N ^ (4.15) 

M=2 
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One can directly read off the (inverse) propagators for the excitation spectrum from this expression 
by evaluating it for p = po- Since f/fe(po) = 0, the — 1 Goldstone bosons are massless and have a 
vanishing decay width as expected. In contrast, the radial mode 0i has a mass mi = ^y2poU'|!{po)/Zl 
and a decay width Fi =71/ ^/2pQU'j!{po). 

4.3 Flow of the effective potential 

Using (4.15) it is straight-forward to derive a flow equation of Uk{p). To that end one simply 
evaluates equation (3.7) for constant argument (p = ^/2p. The result is 



dtUkip)\ 



^U,ip)dtp\^ + dtUkip)\- 
V^pU',ip) + dtU,ip)\^ 



(4.16) 



with 



dtUk{p)\ 



+ 



(AT-l) 



—dtRk- 



(4.17) 



In principle, the frequencies in (4.17) are summed over the discrete, imaginary values po = i^n = 
i2'KTn with J ^ = T as appropriate for the Matsubara formalism. One should not take the 
expression in (4.17) literally at these imaginary frequency, however. In fact the expansion in (4.13) 
was made around a point on the real frequency axis and the function s{po) is only defined there. 
In praxis one performs the Matsubara summation in expressions like (4.17) usually using methods 
based on contour integrals in the complex frequency plane. At the end one has to evaluate residues 
and integrals along branch cuts on the real frequency axis (or close to it for A; > 0) where the 
expansion in (4.13) and the expression in (4.17) can be used. 

In terms of the integral functions defined in appendix A one can write (4.17) as 



dtU{p)\ 



P 2 



U' + 2pU" 7? 7 ^ , 

-'0 I ; 17^1 ^'^i 1: 1^1"' 



fc2 



k 



fU' T 
+ (Ar-l)/of — ,0,l,-,c,rf 



(4.18) 



Up to terms that are proportional to Zp or 7^ (and therefore sub leading), the derivatives of 
(4.18) with respect to p yield 



dtU\p)\ 



„d-2 



{SU'' + 2pU^^^)li('^^^,f,,Z,,lc,d 



+ (Ar-l)?7'7i(^,0,l,pC,rf 



(4.19) 
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and 



dtU"{p)\ 



{3U" + 2pU'^'^) h 



2 ^ fU' + 2pU" -il 



A:2 



T 



. fU' _ T 



(4.20) 



-(iV-l)(f/")'/2(^^,0,l,-,c,f/ 



+ A;2(iV-l) f/(3)/^(^,o,l,pC,rf 



In principle one could attempt a solution of the flow equation (4.17) allowing for a completely 
general form of Uk{p). Numerically one would have to solve a two-dimensional partial differential 
equation for that purpose. Since we are interested in conceptual and qualitative rather than precise 
quantitative results in the present paper, we restrict ourselves to a Taylor expansion of Uk{p) in the 
form ^ 

(4.21) 



Uk{p) = Ukipo) + ml{p - po) + -Xk{p - Po)^ 

the higher orders being neglected for simplicity. In the phase with spontaneous breaking of the 
0{N) symmetry, the (fc-dependent) minimum po,fc is at positive values, po,fc > and the linear 
coefficient ml vanishes by definition. 

In summary, our truncation consists of the /c-dependent renormalized coefficients Zi, jf, and 
Po^k, supplemented by the anomalous dimension rj = —-^dtZ^. 

The flow equation for is obtained by taking the derivative of dtUk{p) in Eq. (4.16) according 

to 

dtXk = g^,dMp)\^^ = 2r^^Xk + ^.f/^p)] _^ (4.22) 

and the last term can be taken from (4.20). In terms of the dimensionless coefficients 7^ = '-/f/k'^, 
Zi, X = Xk, and po = Po,k/k'^ one has 



dtX =2r]iX + 



1 



9A^/2 2Apo,7r,^i,0,c,4 



(iV-l)P/2(0,0,l,0,c,4) 



(4.23) 



The coefficient c is a parameter of the infrared regulator function. A reasonable choice is c ~ 1 and 
the dependence of the final result on the precise value of c can be used for a rough estimation of 
systematical errors connected with the truncation. 

Let us now come to the flow equation for po.fc or po. It is obtained from the condition 



jU'{pQ^k) = dtU'f,{po,k)\p^^ + U'f!{po,k)dtpo,k = 0, 



(4.24) 
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or 

^ ^ (4.25) 

= -V^Po,k - T-9if/'(po,fc)|= • 

In the last equation we have used (4.16). For the dimensionless combination po = Po,fe/^^ this 
results in 



dtpo = - (2 + r]^)po - ^ 



3/i (2Apo,7?,^i,0,c,4) 

(4.26) 

+ (iV-l)/i(0,0,l,0,c,4) ' 



To solve Eqs. (4.23) and (4.26) one needs to specify some initial values for A and p at some large 
value of the infrared cutoff parameter /c = A or, equivalently, at t = ln(fc/A) = 0. One also needs 
the anomalous dimension rj^ and the coefficients Zi and 'jf as a function of t. Equations that govern 
these will be discussed in the subsequent section. 

4.4 Flow of the propagator 

Let us now discuss the (fc-dependent) propagator and its flow equation within our formalism. For a 
constant background field 4>i = 4>o = -\/2po, (t>2 = ■ ■ ■ = 4>n = 0, the propagator is a diagonal matrix 
in momentum space as well as in the space of "internal" or 0{N) degrees of freedom labeled by the 
index M. 

The fc-dependent, renormalized propagator of the radial field 0i is given in the truncation (4.15) 

by 

1 

~ Zi{-pI +p')- is{poH + 2XkPo,k ^^'^^^ 
while the propagator of the Goldstone modes (f)2, ■ ■ ■ ,(pN is 

G2 = (4.28) 

To trace the fc-dependence of these objects one needs besides the flow equations for A^ and po,fc 
also the ones for the discontinuity 7^ and the coefficient Zi as well as the anomalous dimension 
?70 = —Y^dtZ^. To derive these flow equations we expand in the fields according to (4.10). The 
term quadratic in the "fluctuating fields" 50 is in general of the form 

/ \ l^M-Po,-p)Pi{po,p)S(l)i{po,p) + l^(pMi-Po,-P)P2{po,p)S4>M{po,p) \ ■ (4.29) 

"^PO.P [ ^ M=2 ^ J 



fc,2 



Due to the 0(iV) symmetry, the functions Pi{po,p) and P2{po,p) can differ only for p > 0. Within 
the truncation (4.7) their particular form is given by (4.11). We make now an approximation where 
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momentum- and frequency dependent parts proportional to to Zp or 7^ are neglected in the effective 
vertices for the fields 0i and 02 on the right hand side of the flow equation. Using standard methods 
it is then straight forward to derive the flow equations 

1 \ d ^ 
-^dt — ^Pi{qo,0) = 

- [iSpiUy + 2Vf/"f/(3) + 8p3(t/(3))2] A;'^-6 

/go U' + 2pU" U' + 2pU" 7? 71%^^ 
^ \k' P ' P ' fc^' P ' I' 

- 2(iV - V- ^ (f , ^, ^, 0, 0, 1, 1, ^, M 
1 1 9 - 

-4(iv-i)„nV-v(|.^,g.2|.o.z.i,|,c,. 

^9tDisCg(,Pi(go,0) 



(4.30) 



d-6 



= [l8p{Uy + 2Vf/"f/(=^) + 8p3(f/(3))2] A; 

/go + 2pU" U' + 2p^-- 7? 71%^^ 

+ 2(iV - 1) V(f/")^A:^-^ ^ (f ' ^' °' °' ^' ^' P ^' ^) ' 
^9tDisCg(,P2(go,0) 

= 4(iV - 1UU")V-'K (f . ^^:±^, g, |.0,Z.,1, 

We used here the integral functions J and ii" as defined in appendix A (also the operator Disc is 
defined there). From these equations one can easily infer the anomalous dimension 77,^ = —-^dtZ^ 
according to 



= - 4(]V - l)pA' J (0, 2Apo, 0, 7?, 0, Zi, 1, 0, c, 4) . 
In a similar way one can obtain the flow equation for Zi from (4.30) 



(4.31) 



119- 
dtZi =ri^Zi - Pi(go,0) 
Z^ 2go 9go 



i}0=mi 



=770^1 - 18 poA^ J (^go,2Apo,2Apo,7i,7i,Zi,Zi,0,c,4j '^^'^^^ 
+ 2(iV - 1)2 poP J (go, 0, 0, 0, 0, 1, 1, 0, c, 4) . 
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Figure 1. Flow of the interaction strength A 



k- 



According to the definition in (4.14) this should be evaluated as qo = mi/k 
equation for the discontinuity ^\ can be obtained from (4.30) as 



-=—dtDisCqoP2{qo,0)\ 



qo=mi 



+ 18 A^po K (^go, 2Apo, 2Apo, 7f , 7f , ^i, ^i, 0, c, 4 
+ 2(iV - l)2ppo K (go, 0, 0, 0, 0,1,1, 0, c, 4) . 



2Apo/^i- The flow 



(4.33) 



Again this should be evaluated at go = \J'^^Pq/Zi. The second term on the right hand side of 
(4.33), which is due to fluctuations of the massive radial field, will not contribute to the flow 
of the discontinuity since the decay of a massive particle into two particles of the same mass is 
kinematically not possible. However, this is different for the last term in (4.33) which is due to 
fluctuations of the Goldstone modes. A massive particle can decay into two massless ones so that 
a non-zero decay width of the radial field is generated by this term. 

4.5 Numerical results 

Equations (4.23), (4.26), (4.31), (4.32) and (4.33) constitute a closed set of flow equations which 
can be solved numerically. For doing that one also needs to specify some initial values at the UV 
scale = A, of course. For an illustration we choose them here to be 



A(A) = 0.6, po(A) = 0.02, Zi(A) = 1 and 7f(A) = 0. 



(4.34) 



The resulting flow behavior for the 0{N) model in 3 + 1 space-time dimensions is shown in Figs. 1 
to 5. 

In Fig. 1 we show the (rather weak) scale dependence of the interaction strength A. As expected, 
the dependence on k is logarithmic. Strictly speaking, the logarithmic running implies \{k) — )■ 
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Figure 2. Flow of the minimum of the effective potential /3o,fc- 
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Figure 3. Anomalous dimension rj^ as obtained from (4.31). 

for A; — 7- which is cfosely connected to the "triviahty" of the theory. One can clearly separate two 
regimes in Fig. 1. For large scales k the logarithmic running is stronger since both fluctuations of 
the Goldstone modes and of the radial mode contribute. At smaller scales k the contribution of 
the radial mode becomes small since it is suppressed by its non-zero mass mi = A/2Afcpo,fc/^i- The 
transition between the two regimes takes place for mi ^ k. 

In Fig. 2 we show the fc-dependence of the minimum of the effective potential po,fc- For large 
scale parameters it first decreases when k is lowered before it settles to a finite value where it 
remains for — )■ 0. In Fig. 3 we show the anomalous dimension r/^ = —-^dtZ^ as a function of the 
scale. For large k it has an oscillating behavior before it goes to zero for smaller k. 

More interesting for the purpose of the present paper is the flow of the coefficient Zi and the 
discontinuity or decay coefficient 7^ as shown in Fig. 4 and 5, respectively. We show the solution of 
the flow equation evaluated at the frequency go = "^1 as discussed in the preceding section (solid 
lines) but also the corresponding result if the flow equations for Zi and 7i are evaluated at vanishing 
frequency go = instead (dashed lines). Interestingly, Zi changes only very little from Zi = 1 to a 
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Figure 4. Flow of the coefficient Zi as obtained from (4.32) (solid line). We also show the resulting 
behavior if the flow equation is evaluated at go = instead (dashed line). Interestingly, one finds Zi — )■ cx) 
for /c — 7- in the latter case whereas the result is completely regular if the flow equation is evaluated 
on-shell. 
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Figure 5. Flow of the discontinuity coefficient 7^ as obtained from (4.33) (solid line). We also show 
the resulting behavior if the flow equation is evaluated at (70 = instead (dashed line) . As expected, the 
discontinuity 7^ is non-zero on-shell whereas it vanishes for qo = 0. 

slightly larger value if the flow equation is evaluated at go = "^i but diverges for A; — )■ if evaluated 
at go = 0. 

As expected, the flow of 7^ vanishes identically when evaluated at go = but is non-trivial 
at go = mi. One finds that 7^ starts to deviate from zero at some scale /c ~ go, increases then 
strongly before it decreases again slightly and settles to some finite positive value where it remains 
for — )■ 0. The physical origin of this discontinuity in the propagator of the radial mode is that 
it can decay into two massless Goldstone excitations. We emphasize that the present calculation 
takes this decay width in a self-consistent manner into account. 
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5 Conclusions 



We have discussed how functional renormahzation group equations can be analytically continued 
from imaginary Matsubara frequencies to the real frequency axis. For the specific example of a 
relativistic scalar field with 0{N) symmetry we have derived fiow equations for the propagator 
of the radial mode that are evaluated for real frequencies corresponding to the on-shell energy of 
this massive excitation. A prominent and interesting feature is the imaginary discontinuity of the 
inverse propagator which is closely connected to the particle decay width. 

An improved derivative expansion formalism allows to expand in frequencies and momenta 
around their on-shell values in Minkowski space. The thus approximated propagators can be used in 
a self-consistent way in loop expressions on the right hand side of fiow equations if standard methods 
based on complex contour integration are used to perform the Matsubara frequency summations. 
This adapted version of the derivative expansion in Minkowski space is very close to the physical 
dynamics. Since the dynamical properties of an excitation are dominated to a large extend by the 
singular part of its propagator one can expect that the convergence properties are better compared 
to a derivative expansion in Euclidean space. 

We emphasize that our formalism fully conserves Lorentz symmetry. This is facilitated by a 
particular choice of infrared regulator function which decays algebraically both for large spatial 
momenta and large imaginary frequencies. Its rather simple form allows to perform all frequency 
summations or integrations analytically but still leads to expressions for fiow equations that are 
ultraviolet and infrared finite and thus need no further regularization. 

Although we have discussed here for simplicity of notation only relativistic scalar fields, the 
formalism can be extended in a straight-forward way to more complicated theories with bosonic 
and fermionic fields of different spin at arbitrary temperature and chemical potential as well as to 
non-relativistic field theories. 

The formation of bound states or other composite degrees of freedom constitutes a particular 
interesting field of application. Using a recently derived exact fiow equation for composite operator 
fields it is possible to change between a description in terms of fundamental and composite fields in 
a continuous manner during the renormahzation group evolution [19], see also [20, 21]. So far, this 
program was limited to some extend by the difficulty to identify for a given correlation function the 
part that is due to the exchange of some composite degree of freedom. The formalism developed 
in the present paper will allow to determine this part directly from the on-shell singularity for the 
corresponding real frequency. This is in general not possible without the analytic continuation from 
Matsubara space. 

Another interesting prospect for the formalism is the calculation of transport properties. Indeed, 
all quantities that are accessible from linear response theory such as viscosities, conductivities, 
permittivities, permeabilities etc. or the corresponding relaxation times can now be calculated from 
functional renormahzation. The fiow equations of the static and dynamic response functions can 
be determined within a given truncation. In some cases it may be advantageous to introduce 
appropriate composite degrees of freedom explicitly into the formalism. 

In summary, we believe that the analytic continuation of functional renormahzation group 
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equations brings the formalism closer to the physical dynamics in Minkowski space, allows to cal- 
culate more observables and could finally lead to more accurate results with comparatively little 
computational effort. 



A Integral functions 

We define the integral functions according to 



where the frequency po is summed over the discrete imaginary values pq = 2TTiTn with J^^ = T 
T = Tk and t = ln(/c/A). Since the cutoff function is usually proportional to the wavefunction 
renormalization constant Z, the right hand side of (A.l) includes a term proportional tot] = —^dtZ. 
In many situations one has ^ 1 and this correction can therefore often be neglected. 

Note that in terms of a formal cutoff derivative dt which hits only the cutoff Rk, one can write 
for n > 1 

Jpo J \^^) [z{,-pt +P^) + m^k^ - is{pon k + ^Rk\ 
We note also the recursion relation 

d ^ ^ ~ ^ ^ ~ 

7— /„(m^7^z,T,c,rf) = (5on - ri) J„+i(m^ 7^ 2;, T, c, rf). (A. 3) 

For the flow equations of the propagator one needs also the integral functions 

J{qo,m^,m2,'Ji,'y2,zi,Z2,T,c,d) = -k dt 



1 



[zi{-pl+p^) + mjk'^ - is{po)'yfk^ + ^Rk] 

1 

[z2 {-{Po + %kf +p^) + mp2 _ isij,^ + go)7l^^ + ^Rk] 

(A.4) 



and 

d'^-^p 



K{qo, ml, ml, 7i , 72^ ^i^ ^2, c, d) =fc^ '^dt DisCgo 



Po 



{2nY-^ 
1 



X 



[zi{-pl +]P)+ mjk^ - is{po)-ffk'^ + ^Rk] 

1 

[z2 {-{po + qoky mlk^ - is{po + qo)^ik^ + ^Rk] 

(A.5) 
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where the operator Disc is defined (for go > 0) by 

Disc,-J(go) = I [/(go + le) - f{qo - . (A.6) 

For go < the definition of the operator Disc contains an additional minus sign. In praxis the 
integral in (A. 5) has for k > not only a single discontinuity on the real frequency axis but several 
ones on lines that are approximately parallel to the real go axis. They merge only for /c — )• 0. 
Instead of tracing the complete analytic structure, we follow in this paper an approximation scheme 
where the fc-dependent inverse propagator is taken to have only a single cut for real go- The fiow 
equation of this discontinuity is determined by summing the different contributions in Eq. (A. 5). 
This procedure gives the correct result for — )■ and for intermediate k it should be a reasonable 
approximation. 

Using the regulator function given in section 4.1 and the expression for the regularized propa- 
gator given there, one can use standard methods to perform the Matsubara frequency summation 
in Eqs. (A.l), (A. 4) and (A. 5). The integrals along branch cuts get substantially simplified when 
one assumes that they are dominated by the poles that are close to them on the different Riemann 
sheets and can then be done analytically. Some of the remaining integrals over spatial momenta 
can be done analytically as well, others numerically. 
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